Static post-Newtonian limits in non-projectable Horava-Lifshitz gravity with an extra 
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In this paper, we study static post-Newtonian limits in non-projectable Hofava-Lifshitz gravity 
with an extra U(l) symmetry. After obtaining all static spherical solutions in the infrared, we apply 
them to the solar system tests, and obtain the Eddington-Robertson-Schiff parameters in terms of 
the coupling constants of the theory. These parameters are well consistent with observations for the 
physically viable coupling constants. In contrast to the projectable case, this consistence is achieved 
without taking the gauge field and Newtonian prepotential as part of the metric. 
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I. INTRODUCTION 

To quantize gravity in the framework of quantum field 
theory, recently Hofava proposed the Hoi^ava-Lifshitz 
(HL) theory of gravity in which the Arnowitt-Dcscr- 
Misner (ADM) variables [4] are taken as the fundamental 
quantities to describe gravity. By construction, the the- 
ory is power-counting renormalizable, which is realized 
by including high-order spatial derivative operators [up 
to six in (3+l)-dimensional spacetimes]. The exclusion of 
high-order time derivative operators, on the other hand, 
ensures that the theory is unitary, a problem that has 
been faced in high-order derivative theories of gravity for 
a long time Q . Clearly, this inevitably breaks the general 
diffeomorphisms , 



feM^_C^(t,a;),(M = 0,l,2,3). 



(1.1) 



Although such a breaking in the gravitational sector is 
much less restricted by experiments/observations than 
that in the matter sector 0, 01 j it is still a challenging 
question how to prevent the propagation of the Lorentz 
violations into the Standard Model of particle physics Q . 
Hofava assumed that such a breaking only happens in 
the ultraviolet (UV) and down to the foliation-preserving 
diffeomorphism. 



St 



-fit), Sx' = ~C(t,x),{t^ 1,2,3), 



(1.2) 



often denoted by DifF(M, J^). In the infrared (IR), the 
low derivative operators take over, presumably providing 
a healthy low energy limit. 

The breaking of the general diffeomorphisms immedi- 
ately results in the appearance of spin-0 gravitons in the 
theory, in addition to the spin- 2 ones, found in general 
relativity (GR). This is potentially dangerous, and leads 



to several problems, including instability, strong coupling 
and different speeds of (massless) particles To resolve 
these problems, various models have been proposed 
including the healthy extension of the non-projectable 
HL theory ^8|. In the healthy extension, the instability 
problem was fixed by the inclusion of the term /Soflja* in 
the gravitational action, where /3o must be in the range 
< /3o < 2, and a.^ = N^i/N with TV being the lapse func- 
tion in the ADM decompositions . The strong coupling 
problem is resolved by introducing a new energy scale 
M*, so that < = y^Mpi, where M* denotes the 
suppression energy of the high-order spatial operators, 
A„ the would-be strong coupling energy scale, and Mpi 
the Planck mass. Clearly, in order for this mechanism to 
work, one must assume that /3o 7^ 0. Observational con- 
straints requires 10^° GeV < < lO^^ GeV [1]. The 
low bound was obtained by assuming that also sets 
the suppression energy scale in matter fields, while the 
up bound was obtained from the preferred frame effects 
[3, [l3| • In addition, to avoid the Cherenkov radiation, 
one must require that the speed of the spin-0 gravitons 
be superluminal [llj ^. 

A more dramatical modification was proposed recently 
by Hofava and Melby-Thompson (HMT) [i3|, in which 
an extra local U(l) symmetry was introduced, so that 
the symmetry of the theory was enlarged to. 



[/(I) X Diff(M, J"). 



(1.3) 



Such an extra symmetry is realized by introducing a 
gauge field A and a Newtonian prepotential tp. Because 
of this extra symmetry, the spin-0 gravitons are elimi- 
nated [l3, [l3|- As a result, all the problems related to 
them, such as the instability, strong coupling and dif- 
ferent speeds in the gravitational sector, as mentioned 



^ This can be easily seen using the equivalence between the 
'Electronic address: llk3141 59@hotm ail.coml Einstein-aether theory and the nonprojectable HL theory in the 

t Electronic address: | anzhong' wang @baylor.edu| IR limit d [3 • 
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above, are automatically resolved. This was initially done 
with the projectability condition N = N{t) and A = 1 
, where A characterizes the IR deviation of the theory 
from GR. It was soon generalized to the case with any 
value of A ilSil . in which the spin-0 gravitons are still elim- 
inated |l5l.ll6l| . Although the strong coupling problem in 
the gravitational sector disappears, it still exists in the 
matter sector [16], but can be resolved also by introduc- 
ing a new energy scale M*, so that M* < |A — 1| ^/*Mp/ 
[17| . Cosmological applications of this model were con- 
sidered in [Tsl. Il9l|. and found that it is consistent with 
current observations [20]. On the other hand, the stud- 
ies of solar system tests recently showed that the theory 
is consistent with observations if and only if the gauge 
field A and the Newtonian prepotential tp are part of the 
metric [U 2. 

A non-trivial generalization of the enlarged symmetry 
()1.3p to the nonprojectablc case N — N{t, x) was re- 
cently worked out in [26,, 21} , and showed that the only 
degree of freedom of the model in the gravitational sec- 
tor is the spin- 2 massless gravitons, the same as that in 
GR. Because of the elimination of the spin-0 gravitons, 
the physically viable region of the coupling constants is 
considerably enlarged, in comparison with the healthy 
extension [8], where the extra U(l) symmetry is absent. 
In particular, the requirement /3o 7^ now is dropped out. 
Furthermore, the number of independent coupling con- 
stants is also dramatically reduced, from about 100 down 
to 15. The consistence of the model with cosmology was 
showed recently in [27l - [29| . and various remarkable fea- 
tures were found. 

In this paper, we study the consistence of the model 
proposed in [2^, [l^] with the solar system tests. In par- 
ticular, we investigate its static post-Newtonian limits. 
The paper is organize as follows: In the next section 
(Sect. II), we give a brief introduction to the model, 
and in Sec. Ill we find the spherical static solutions in 
the IR limit, from which we can see that the Birkhoff 
theorem is not applicable to the current theory, because 
of the breaking of the Lorentz symmetry. As a results, 
there are several classes of static asymptotically fiat vac- 
uum solutions. Then, in Sec. IV, we study the post- 
Newtonian limit for each class of these solutions, and find 
the Eddington-Robertson-Schiff parameters, 7 and /?, in 
terms of the coupling constants of the theory, where 7 is 
related to the amount of spatial curvature generated by 
the spherical source, and f3 to the degree of non-linearity 
in the gravitational field. By doing so, we show that the 
observational constraints on 7 and /3 can be easily satis- 
fied within the physically viable region of the phase space 
of the coupling constants of the theory. This is true with- 



^ When A = 1, the theory is consistent with the solar system tests 
even without the gauge field and the Newtonian prepotential be- 
ing part of the metric |22| , |23( , but the corresponding cosmology is 
quite different from the standard one [16| |. Other considerations 
of this model can be found in 1241. |25|| . 



out taking the gauge field A and Newtonian prepotential 
if as part of the metric, in contrast to the projectable 
case [21|. In Sec. V, we summarize our main results and 
present some discussions and remarks. 



II. NON-PROJECTABLE HL THEORY WITH 
U(l) SYMMETRY 

In this section, we shall give a brief review of the non- 
projectable HL theory with the enlarged symmetry (|1.3p 
p6l \2T\ . The fundamental variables of the theory are 
(TV, N^, gij, A, ip^, where N'^ and gij are, respectively, 
the shift vector, and 3-metric of the leaves t = Constant 
in the ADM decompositions Under the local U{1) 
symmetry, the fields transform as 



6aA — a — N^Via, S^ip — —a, 



(2.1) 



where a is the generator of the local U{1) gauge sym- 
metry, a = da/dt, and Vi the covariant derivative with 
respect to the 3-metric gij. Under the Diff(Af, T), they 
transform as, 

6N = (''VkN + Nf + Nf, 

Sgij = VjCj + + fgij, 
6 A = C^.A + fA + fA, 

dip = f^ + Cy^V■ (2.2) 

The general action reads [1^ [l^j , 



s = e 



dtd 



Ck — C.\ 



La + Lu 



M 



(2.3) 



where = l/(167rG) with G being the Newtonian con- 
stant. Cm describes matter fields, and 



Cv 

Ca 
r 



C^, 



= 4(2A,-i?), 
= ipQ"-^ {2Kij + S/iS/jip + ai\/jLp) 
+ (1- A)[(Av7 + a,VV)' 

+ 2(Aip + aiV'-ip)K 
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4{V,Vjip) a(^tSJi)Lp 
+ 5 {a(iVj)ip) a(^k^^Lp 
+2{\/i^.,(p)aj)i^k^i)ip + 6Kijai^iWk)(p , (2.4) 



with A 



V^, and 
1 



2TV 



i~g^j+V^Nj+VjN,). 
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72^^ + 73%-R'" , 75 



/•a 



1 



/3i {a,af +(32{a\f 



+/33 (ajfl*) a-'j- + /34a'-' a„ + /Jg (0^0') i? 

1 



+ -/?8(AaO . (2.5) 



Here i?^^-, i? are, respectively the Ricci tensor and scalar 
of Qij, and Cij denotes the Cotton tensor, defined by 



V5 ^ 4 ' 



(2.6) 



with e^^'^ = 1, etc. A,7„,/3s and Kg are the coupling 
constants of the theory. In terms of Rij, we have |27| . 



2 2 ^ 8 



(V,i?,fc) V'i?^" +VfeG 



where 



2 8 
The infrared limit requires that 

A=i7oC', 71 = -1, 



(2.7) 



(2.8) 



(2.9) 



where A denotes the cosmological constant. 

The variations of the action with respect to N, N^,A,ip 
and gij, yield, respectively, the Hamiltonian, momentum, 
A-, ip- constraints, and dynamical equations ^, given ex- 
phcitly in [13 . 



III. SPHERICAL STATIC VACUUM 
SOLUTIONS 



In this paper, we consider the s phe rically symmetric 
static spacetimes, described by [23l l30j|. 

N^N{r), N'^h{r)5i, A = A{r), 



W) 



+ r^dn^, Lp^Lp{r), (3.1) 



^ In the expression of Fa^ given by Eq.(3.22) in [27| . the coefficients 
/3s should be replaced by /9s, where /3s = (— /3o,/3n) with n = 
1, 2, 8, while 79 defined in Eq.(3.24) should be defined as 79 = 
75, instead of 75/2. 



where dfl^ = dO"^ +sin^ Odcf)^ . Using the U(l) gauge free- 
dom, without loss of the generality, we set 



= 0. 



(3.2) 



Then, we find that C^p = Q. 



In the IR, the spatial curvature is very small, and its 
high-order derivatives are negligible, so we can safely as- 
sume that Cv has only three non- vanishing terms, given 

by 



Cv = 2A - R - ^3oa^a\ 



(3.3) 



In the solar system, the effects of cosmological constant 
are negligible. In addition, the curvature is extremely 
small. So, without loss of generality, we can safely set 
A = Ag = 0. Then, we find that in the present case there 
are only four-independent equations, which can be taken 
as the Hamiltonian and momentum constraints, the A- 
constraint, obtained from the variation of the gauge field 
A, and the rr-component of the dynamical equations, 
given, respectively, by 



/ h J ^{f h rj ^ 




f 2N' 

■2 
•2 



f r) N 



h' 



0, (3.4) 



/" /'^ f'N' A fN' 1 



(rf)' = 1, 
A 



fl/2 



2hfr_2N^\ 
r \ f ' 



Qir) 




where 



Qir) = 



+ 



+ 



f \ f N 

Ah' fh' 2N' 
~ \ h^~ 

N f2N' f 



/ h 



16 

r 



2/1/2 



EL 
1 

" 2^ 



fj 4/1/27V 



N'^. (3.8) 



Therefore, in the present case, we have four independent 
ordinary differential equations for four unknowns, /, h 
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and A. In particular, Ea. p.6|) has the general solution, 

(3.9) 



where B is an integration constant. Once / is known, 
from Eqs. (l3.4p and (13. 5p one can find TV and h, while 
Eq.(E31) yields. 



2B 



Ao + / Q{r)dr , (3.10) 



where Aq is an integration constant. Therefore, the main 
task now reduces to solve Eqs. p.4p and p.5|) for N and 
h. In the following, let us consider the two cases h = 
and h ^ 0, separately. 



h = 



When f3o ^ 0, from Ea. (l3Tl|) we find that 

\ 2 

2B\ 



N{r) = Na + Ni\ 1 



(3.14) 



where A^i is another integration constant. Inserting it 
into Eqs. (|3^ and ([330]) . we obtain 



A{r) = n'^+A,\Ii-^-N!^ 

B 



-iPo - l)Nt 1 



B. /i / 



(3.15) 



In this case, to solve Eqs. (l3.4p and p.Sp . we consider 
the cases A = 1 and A 7^ 1, separately. 



In this case, the momentum constraint is satisfied iden- 
tically, while Eq. p.4p reduces to 



f3o 



EL 



4\ N' 
r) 77 



0. 



(3.11) 



When f3o =0, it is satisfied identically, while Ea. (|3.10p 
yields. 



N^Nox 1 



2B 

r 



Air), 



(3.12) 



where A'o is a constant, and A{r) is undetermined. In 
particular, when A = 0, the above solutions reduce to 
the Schwarzschild solution. 



/ 



7V2 



1 



— , A{r) 



0. 



(3.13) 



Note that in writing the above expression, we had set 
A^o = 1 by rescaling t. It is remarkable to note that 
the above solutions are independent of the coupling con- 
stant A. That is, in contrast to all other models of the 
HL theory Q, the Schwarzschild solution (I3.12p is a so- 
lution of the HL theory not only for A = 1 but also for 
any value of A in the IR limit. In addition, in the non- 
projectable HL theory without the U(l) symmetry, to 
solve the instability problem of the spin-0 gravitons, /3o 
is necessarily non-zero, /3o 7^ 0. Otherwise, the spin-0 
gravitons become unstable. However, in the current case 
the U(l) symmetry eliminates the spin-0 gravitons [27| . 
so the instability problem is out of question, and the 
physically viable region of the phase space includes the 
point f3o — 0. Furthermore, in the projectable case, we 
have N = N{t), and the solutions with h — take the 
Schwarzschild form only when the gauge field A and the 
Newtonian pre-potential ip are part of the metric [l^ [2l| . 
The above observations are important, when we consider 
the solar system tests. 



1. A = 1 

Depending on whether TV is a constant or not, we can 
have different solutions. In particular, when N — Con- 
stant, say, iVp, without loss of generality, we can always 
set it to one by the rescaling, t — t/Nf). Then, from 
Eqs.dXil) and ([311), we find f'h = 0. Since /i 7^ 0, we 
must have /' = 0, or S = 0. Then, Eqs. (|Ot yields, 
(r/i^) = 0, which has the general solution. 



(3.16) 



where is an integration constant. This is precisely 
the Schwarzschild solution, but written in the Painleve- 
GuUstrand coordinates [3l|. Inserting it into Eas. (l3.8p 
and p.lOp . and considering the fact that now we have 
N = f = 1, we obtain 



A(r) - Ao, 



(3.17) 



where Aq is a constant. 

When N is different from a constant, Eg. p.Sp has the 
general solution. 



Then, Eq.daH) and Eg. ({310)) yield. 



Hr) = No 
A{r) = — 



r-2B 
12r2 



Ci + 3BPo In 



r-2B 



(3.18) 



1/2 



8r./l-¥ 



2B - [r - 2B) In 



-Ao\\ 



2B 



r-2B, 

(3.19) 



where Ci is a constant. 
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2. A / 1 

In this case, let us consider the two cases, B — Q and 
S 7^ 0, separately. 

Case (i) B = 0: Now Eqs. (03]) and reduce, 
respectively, to 



ehl' 

^2 









A* r 


\^ 




r 


7?" 


f 4(1 




+ 2(l + 2e) 



= 0, 
(3.20) 



where e = A — 1 . 

When /3o = 0, the solutions is given by 

h{r) = Hor°'±, N{r)^Nar^±, 
Air) =. Ao+Ei{r)+eE2ir), 



(3.21) 



(3.22) 



where Aq and Hq are constant, and 

a± = -2- 



2 ± V4 + 6£ 



/3± 
Ei{r) 



-2 + a±+ a\ 
2 + (2 + a±)e^ 

4A^or^.(2.,-/3,) k°-^'^(^"^'^^) 
X (4 + /3±/3o) - 2i/2r2"± (1 + 2a± + 2/?±) 



i?2(r) = 



2 (2 + a±)(Q;± 



2/?±) 



4iVo(2a± - /3±) 



,2a±-,9± 



(3.23) 



It is interesting to note that, taking the "-" sign in the 
above solutions, they reduce to the Schwarzschild solu- 
tion when e — > 0. This is in contrast to the case with the 
projectability condition N = N{t) [12], in which it was 
found that such relativistic limit does not exist. 

When /3o 7^ 0, it is found difficult to obtain exact so- 
lutions. Instead, we consider the cases where e is very 
small, as one would expect that physical viable solutions 
must be very close to that of GR where Xor = 1- Thus, 
expanding h{r), N(r) and A{r) in terms of e, we find that 



h{r) 

N{r) 

Mr) 
Hence, 



± 



Hp ^ 9{N^l3or + rgln{r)) 



+0(e^), 



iVi - -iVo ln(r) 



9 



Ao+e--7Vo(/3o-2)ln(r) + 0(e2). 
16 



_9_ 

16' 



(3.24) 



(3.25) 



as r — 7^ cxD. Thus, unless /Jqc = 0, the solutions are 
not asymptotically flat. When /3o = 0, it reduces to 
the last case, and the corresponding solutions are given 
by Eq. (|3.22p . When e = 0, the corresponding solutions 
are given by Eas. (|3.18p and (|3.19l) . Therefore, the case 
e/?o 7^ has no physically viable solutions, and in the 
following we shall not consider it any further. 

Case (ii) i? 7^ 0: In this case, to obtain exact solutions 
is found also very difficult, and instead we expand h(r\ 
N(r) and A(r) in terms of e, and find that 

N[r) = iV(r) +eiV(r) +C'(e2), 
h{r) = h{r)+eh{r)+0[e^), 

A{r) = i(r)+eV7 + (3-26) 



where /i(r), iV'(r) and A{r) are the solution for e = 0, 
which are given by Egs. p.lSp and (|3.19p . while 



f 



h{r) = 



rh{r) 



Hi 



ha{r)di 



{3B-r)h' 
2h " {2B - r)h 

- 8Br + 2r2 



2r(r - 25)2 ' 

2 

hair) = Npppr^ {^^N" + N') - ^{hy 



Nofi 16/3 
No{l-5f) + 8^N 
fioNpB^N 



4iVo/2 



Qir) 



2r2/l 

No{5 - 2/ - 35/2) ^ _ 2rN')fi 



1 



, , —hh-hh' + -h"h 

Nop V r f 2 



rf 

B N 



2Norp N^fi 
N' h' 

l3oB 



h'h 



BN 
r/2 



^ " Yo/2 V A J 
^— (BN + 2rJN') . 



(3.27) 



From the above expressions, we can see that the solutions 
are also not asymptotically flat. In particular, we have 
N{r) cx ln(r). Therefore, this class of solutions is also 
physically discarded. 
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IV. SOLAR SYSTEM TESTS 

In the solar system tests, the metric is usually cast in 
the diagonal form, 



(4.1) 



where 5* and $ of the gravitational field, produced by 
a point-like and motion- less particle with mass M, are 
parametrized in the forms. 



^2* 



„2* 



= l-2x + 2(/3-7)x'-f 
= 1 + 27X+..., 



(4.2) 



where /3 and 7 are the Eddington parameters. For the 
Schwarzschild solution, we have Pgr = Jgr = 1- 

In the solar system, we have rg = GMq/c^ ~ 1.5 km, 
and its radius is tq ~ 1.392 x 10^ km. So, within the 
solar system the dimensionless quantity x[= GM/{rc'^)] 
in most cases is much less than one, x ^ '^g/'^O — 10"^- 
The Shapiro delay of the Cassini probe [sl] , and the solar 
system ephemerides [11] yield, respectively, the bounds 



7-1 = (2.1 ±2.3) X 10"^ 
/?-l = (-4.1 ±7.8) X 10"' 



(4.3) 



To fit the solutions found in the last section with above 
observations, let us consider the cases /i = and h ^ 0, 
separately. 



A. h = 

In this case, let us first consider the solutions given by 
Eas. (|3J)) and (|3.12p . Expanding A{r) in the form. 



A{r)^A,x + A2x' + 0{x') 



(4.4) 



we find that 



1 



N'ir) = l-2{c7-Ai)x + 2l^A2~c7Ai + -Aijx' 

+o{x'), 

r\r) = l + 2ax + 4a\2 + 0(x'), (4.5) 

where a = ^jj- Note that in writing the above expres- 
sions, we had set A^o = 1 by rescaling t and A. Comparing 
Ea. (|i3)) with Eq.dl^I), we find that 

Ai = 7 - 1 < 10"^ 

A, = (/? - 7) + il±ll (^ _ 1) < 10-5, 

a = 7 ~ 1 + (2.1 ± 2.3) X 10"^ (4.6) 

For the solutions given by Eqs. (l3.9p and (|3.14p . we have 

N\r) = {No + Nif-4<jNi{No + Nifx 
+2a^Ni {2Ni ~ No) {Nq + Nif x^ 
+0 (x') , (4.7) 



where /(r) is still given by Ea. (l4.5p . Then, comparing 
them with Eq. (|4.2l) . we obtain 



2/3-7=2' 

' 27 ' 27 



(4.8) 



Clearly, this class of solutions does not satisfy the solar 
system tests, and must be discarded. 



B. /i / 

To study the solar system tests for this class of solu- 
tions, we need first transform the experimental results of 
Eas. (|4?2)) and (|43)) into the form. 



ds^ = -N'\r)dt^ + {dr + h(r)dt] ^ + r'^dfl^, (4.9) 
which can be done by the coordinate transformations, 

r = f-y" e-*v/e2* - /-^dr, (4.10) 
where / is given by Eq. p.9p . and 



iV = e*+*v^, /i = e*V/2e2*-/, (4.11) 
or inversely. 



„2* 



7,2 \ -1 A^2 

^*-'A^^-yj (4.12) 



It should be noted that, although the coordinate transfor- 
mations (j4.10p are forbidden by the foliation-preserving 
diffeomorphisms, we assume that experimental results 
can be expressed freely in any coordinate systems. 

Now let us first consider the solutions given by 
Eas. ((XTB|) and ([XT?)) with 1 ^ /. Clearly, this 

is the Schwarzschild solution. 



2* -2* 1 o 

e = e = 1 - 2x, 



(4.13) 



with /3 = 7 = 1, as it is expected. 

On the other hand, for the solutions given by 
Egs. iprS)) and ([XT^ . we find that 



e^* = l-(2+^)ax-/3oyX^ + 0(x^), 



„2* 



l+(2+|^)-X + 0(x^) 



(4.14) 



Note that in writing the above expressions, we set A^o = 
1. Comparing the above expressions with Eq. (|4.2p . we 
find that 

(GM 

7=1, Ci = 24f — -B 



-i=-r 



V 2GAf 



(4.15) 
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Clearly, by properly choosing /3q and B, the conditions 
(|4.3p can be easily satisfied for this class of solutions. 
For the solutions p.22|) . we find that 



„2* 



B \ 



20± 



Hi 



B 



„2* 



B 



2/3± 



B 
<^X 



2/3± 



(1 - 2ax) 



(4.16) 



Comparing the above with Eq. 



1 



we find that 



(4.17) 



This means A = 1. Thus, in this case the solutions are 
consistent with the solar system tests only when A = 1. 



V. CONCLUSIONS 

In this paper, we have studied the observational con- 
straints of the nonprojectable HL theory with the en- 
larged symmetry (jl.3p within the solar system. In partic- 
ular, we have first found the static spherical vacuum so- 
lutions of the theory in the IR limit, and then considered 
their post-Newtonian limits. Because of the breaking of 
the Lorentz symmetry of the theory, the BirkhofF theo- 
rem does not hold here, and there exist several families of 
vacuum solutions that are asymptotically flat and satisfy 
the observational constraints imposed on the Eddington- 
Robertson-Schiff parameters, 7 and /3, given by Ea. (l4.3p 
in the spherically symmetric spacetimes. In contrast to 
the projectable HL theory with the enlarged symmetry 



(II. SP jl3l-ll6l|. the consistence is achieved without the 
gauge field and Newtonian prepotential being part of the 
metric 

A remarkable feature is that the Schwarzschild solution 
p.l3|) written in the Schwarzschild coordinates is also a 
physically viable solution of the nonprojectable HL the- 
ory with the enlarged symmetry (11.31) in the IR not only 
for A = 1 but also for any value of A. This is different 
from all other models of the HL theory proposed so far 

a. 

Applying the general results presented in |3S] to the 
physically viable spherical vacuum solutions obtained in 
this paper, one can easily construct slowly rotating vac- 
uum solutions of the nonprojectable HL theory with the 
enlarged symmetry (jl.3|) . Some of such constructed solu- 
tions clearly represent slowly rotating black holes in the 
IR limit [li]. 

Post-Newtonian approximations are generally charac- 
terized by ten parameters Q. The studies of observa- 
tional constraints in the nonprojectable HL theory with- 
out the local U(l) symmetry [8'| showed that the preferred 
frame effects imposed the most stringent constraints on 
the suppression energy M*. The structures of the non- 
projectable HL theory with and without the local U(l) 
symmetry are quite different. In particular, the spin- 
gravitons do not exist in the current model. Thus, it 
would be very interesting to see which kind of constraints 
the preferred frame effects can impose on as well as 
on other coupling constants of the theory. We wish to 
come back to this issue soon in another occasion. 
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